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I. INTRODUCTION 


The Atmospheric Boundary Layer (ABL) is defined by Stewart (1979) as that 
portion of the lower atmosphere which has turbulent flow and is in direct contact with 
the earth's surface. The ABL extends from the carthisesuriace to a helentmol sever 
ten’s to thousands of meters. All momentum, water vapor, and thermal energy 
exchanges between the atmosphere and the earth’s surface take place within the ABL. 
The stability of the ABL is dependent upon its density structure. The density structure 
in turn 1s a function of water vapor distribution and whether the layer is heated or 
cooled from below. While both stable and unstable ABLs have some energy for 
vertical mixing as a result of shear production of turbulent Kinetic energy (TKE), when 
there is heating from the earth’s surface, the heating produces a density distribution 
which enhances vertical mixing. In this unstable case, the overall level of TKE, and 
hence vertical mixing, 1s greater. The unstably statified boundary layer laver can be 
further subdivided into three separate layers: the surface layer, the mixed layer, and a 
capping inversion layer. 

For the near-neutral to unstable case, depending upon the magnitude and 
direction of the thermal energy flux, the depth of the surface layer typically ranges 
from a few meters to several ten’s of meters, and is conventionally considered to 
occupy the lower tenth of the ABL (Wyngaard, 1973). Phvsical continuity between the 
atmosphere and the earth dictates that wind velocitv be zero at the immediate surface 
over land and equal to the surface drift current over water. Above the surface, velocitv 
Increases in approximately a logarithmic fashion to the top of the layer. Thus, strong 
vertical wind shear 1s common in this layer. It, along with convection, produces 
turbulence which facilitates the transport of moisture and momentum within the layer. 
The relatively homogeneous mixed layer directly above the turbulent surface layer 
constitutes the majority of the unstable ABL. Turbulence tends to continuously erode 
local gradients and, thus, conditions are well-mixed up to the inversion layer. The last 
50 to 100 meters of the convective ABL is usually an inversion layer. In it the kinetic 
energy available for mixing is damped by strong stable vertical temperature gradients. 
The gradients have a “capping” effect on the convection below and effectively suppress 


exchange of physical properties with the free atmosphere above. 
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If the ABL is cooled from below it will be stable. Turbulence is suppressed by 
the resulting density structure, and generally little muxing will occur except in the 
surface layer. All exchanges between the atmosphere and the earth’s surface are thus 
greativeredwecd. 

It 1s in the marine atmospheric boundary layer (MABL) that the oceans and the 
atmosphere exchange energy directly in the form of turbulent heat, moisture and 
momentum fluxes. Gases and sohds are also in continuous exchange between cach 
fluid. Wave activity can toss small water droplets directly into the lower atmosphere 
where they become suspended and eventually evaporate. Normally, evaporation from 
the sea surface creates a shallow layer of rapidly decreasing relative humidity from 
100% at the sea surface to a lower value (80-90%) directly above the sea surface. This 
low level feature produces an extended radar propagation phenomenon known as an 
evaporation duct because of its formation mechanism. The evaporation duct extends 
only tens of meters above the ocean’s surface but is of great importance to the mariner 
because it occurs in varying degrees throughout all oceanic areas and results in over- 
the-horizon radar propagation. 

The ABL over land shows much more variation than the MABL. Seasonal and 
diurnal changes are both more pronounced. Larger variations in the vertical 
temperature structure occur over land than over the ocean but there is usually less 
vertical humidity variation. Radiative heating and cooling can cause large diurnal 
fluctuations in the height of the ABL from near surface values at might to over 2000 
meters during the dav. 

Propagating electromagnetic (EM) waves, unless in a completely homogeneous 
medium, will experience some degree of bending due to changes in the index of 
refraction. The earth's atmosphere is normally a very inhomogeneous fluid. Certain 
regions such as the ABL characteristically have large mean gradients in temperature 
and or humidity. Rapid vertical changes in both temperature and humidity create 
lavers which significantly refract propagating EM signals. This phenomenon is readily 
apparent, for example, in the evaporation duct at the base of the MABL and in the 
inversion layer at the top of the ABL. Refraction of EM waves will receive more 
cetailed discussion in Chapter two. 

The Integrated Refractive Effects Prediction System (IREPS) is a shipboard 
environmental data processing and display system which is used to predict the effects 


of refraction on electromagnetic signals for naval surveillance, communications, 
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electronic warfare, and weapons guidance systems. Environmental data usuallv 
cbtained from radiosondes are processed on a desk top computer, the Hewlett-Packard 
S845, to derive a comprehensive assessment of refractive effects in the lower 
atmosphere. Knowledge of the refractive environment can be used to maximize tactical 
advantage. The most recent IREPS, version 2:2yisworesently sscheduledmtcmbemam 
integral part of the Tactical Environmental Support System (TESS)(NOSC TD 659, 
1981), soon to be introduced into the flcet. 

The performance of the [REPS micdel 1s “dipectily velateds tom nemenialin mania 
timeliness of atmospheric information entered into the program. Accurate and timely 
data must be used if the resulting coverage displays are to indicate the true refractive 
structure of the atmosphere. Knowledge of the immediate environment and its 
dynamical evolution are necessary for a proper interpretation of IREPS results. The 
IREPS revision 2.2 User's Manual on page 49 (NOSC TD 659, 1981) states that the 
assumption of a horizontally homogeneous atmosphere is valid about 85 percent of the 
time for the purpose of making refractive effects assessments. Given the natural 
variability evident in atmospheric soundings made with measurement systems having 
high temporal resolution (e.g. Webb, 1977), this statement may be too optimistic. 
Thus, one of the objectives of this thesis is to examine the validity of the assumption of 
the horizontally homogeneous atmosphere. By selecting a data set in an oceanic area 
which was previously believed to be homogeneous, we investigate how the vertical 
structure of the inmimediate atmosphere controls the nature of radar propagation. We 
are concerned not only with the spatial variability of the atmosphere, but also its 
temporal variability. The atmosphere undergoes continous dynanuc evolution in 
diurnal as well as synoptic time scales. The IREPS User’s Manual notes that “the user 
Should be aware of the changing state of the atmosphere and try to acquire and use 
refractivity measurements that are appropriate to the planned time and place of 
pertinent operations.” Our second objective is to attempt to qualitatively deterniine the 
length of time an IREPS product could remain an operationally useful predictor of the 
refractive structure of the surrounding environment. Thirdly, we expect that in order 
to derive an accurate IREPS assessment of the refractive structure of the atmosphere, 
high vertical resolution data are required. Presumably, unless the vertical resolution of 
the input data is sufficient to define relatively thin lavers of strong temperature and 
humudity gradients which could affect atmospheric refractivity (e.g. the inversion layer 


of an ABL), the resulting IREPS output will not be an accurate representation of the 
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actual atmospheric refractive structure. Accordingly, the differences in IREPS output 
products, depending solely upon the quality of resolution of the input data, will be 
compared. 

The FORTRAN code used in TESS was used to run various [REPS products at 
tne Naval Environmenta! Research Prediction Facility (NEPRF). Although TESS ts 
not Vet operational, the newest revision of the IREPS code, version 2.2, was used for 
all the atmospheric refraction computations. Thus, the results of this work should be 
directly applicable to any naval unit depending upon TESS for atmospheric structure 


and refractive effects information. 
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It. ATMOSPHERIC REFRACTIVITY 


The index of refraction of a medium) n. 1s defined By (2,1) where c 1s ene velocms., 


of an EM wave in free space and v 1s the velocity of the same EM wave in the medium. 


n= Cy (2.1) 


For the atmosphere, the refractive index is defined to be the ratio of the velocity 
of propagation of an EM wave in a free space to that in the air. Since the velocity of 
an EM wave in a free space is always faster than that in any medium, the index of 
refraction 1s alwavs greater than one. In the standard atmosphere, a typical near sea 
level value for n 1s approximately 1.0003 (Bean and Dutton, 1966). For convenience, 
the refractivity N is defined (2.2) such that normal atmospheric values of N units 
range from 250 to 400. 


N = (n-1)* 10° (2.2) 


As Stated by Battan (1973), “In dry air the index of refraction has the same value 
over almost the entire range of the electromagnetic spectrum: it 1s the same for light 
and radio waves. However, when water vapor is added to the air, the value of N for 
the mixture becomes frequency dependent. It is well known the water molecule 1s 
polar in nature and that the dipole moment of the molecule has a different response to 
different-frequencyv radio waves. With the extremely high frequencies of visible light, 
the water molecules are electronically polarized. With the lower frequencies of radar 
waves, the water molecules not only acquire electronic polarization, but also reorient 
themselves rapidly enough to follow the electric-field changes. As a result, the index of 
refraction (and dielectric constant) of water vapor is greater for radio than for optical 
frequencies.” Temperature, water vapor and pressure are the major variables of the 
atmosphere which determine its refractivity. Bean and Dutton (1966) express 
refractivity in terms of temperature T in Kelvins, with pressure P and water vapor 


pressure e in mb as follows in (2.3): 


N = 77.6 P/T - 5.6 e/T + 3.75 * 10° e/T? (2.3) 


14 


In an atmosphere of constant N units, no bending of an EM wave could occur 
regardless of the value of N. EM refraction is dependent upon the gradients of N. 
Since gradients of pressure, temperature and humidity occur throughout the 
atmosphere, it follows that gradients of N must also exist. It can be shown (Battan, 


Pee iewmicn tie Sradicnt Of Vie. d.\/dz) is equal to -157 Km} 


, a propagating 
EM wave will bend with a curvature exactly equal to that of the earth. This would 
cause a horizontally propagating EM wave to remain constantly parallel to the earth’s 


I would 


surface always at the same height. Any value of dN/dz less than -157 Km” 
cause an EM wave to bend with greater curvature than the earth’s surface. Therefore, 
-157 Km’! is the threshold for trapping of an EM wave. 

It 1s sometimes convenient to portray the earth’s surface as being flat and to 
present the refraction of an EM wave in this frame of reference. This can be done 
simplv by mathematically subtracting the earth’s curvature from the EM wave and the 
earth. A modified index of refraction, M, 1s derived on this basis. It is often used in 
vertical profiles to provide a direct indication of refracting conditions that would cause 
= le a CeO bemdecitiier towardeor way trom the earths surface. (2.4) (after 
Battan, 1973) relates the modified index of refraction at any height to the refractivity at 


that height. 


M = N + 157z (2.4) 


In vertical displavs of the modified refractivity, any layer where M decreases with 
increasing height (a negative value of dM/dz) is an indication of a region where an EM 
wave would bend with a curvature greater than that of the earth. Depending on the 
height, strength and extent of this type of layer, an EM wave could be refracted into 
the earth’s surface. Subsequent reflections would direct this EM wave back into the 
atmosphere such that over-the-horizon radar ranges can result. This is the 
characteristic of an evaporation duct which makes it so important in naval operations. 
Elevated trapping layers, too, can result in extended radar ranges. However, depending 
upon the height of the particular layer, this form of ducting may be of more 
importance to airborne radar svstem than to surface based units. Fig. 2.1 portrays the 
basic types of ducting that may occur due to typical M vertical profiles. 

Profile (a) of Fig. 2.1 depicts an elevated duct which is the type often found in 


an inversion layer. This type of ducting condition can occur when warm, dry air exists 


over a layer of cooler, moist air. This phenomenon 1s particulary pronounced over 
regions in eastern oceans and to a lesser degree in regions with trade wind inversions. 
The top of the duct is defined as the height where M reaches its minimum value. The 
thickness of the duct is found by dropping a vertical line from the top of the duct down 
towards the surface until the M profile is again intersected. Duct strength is defined as 
the maximum range of M values within the limits of the duct. The optimum coupling 
height (CCH) is the height where the dM/dZ profile changes from a positive to a 
negative value. Profile (b) is an example of a surfaced-based duct. These ducts can be 
formed by warm, dry air being advected over a cool body of water, or by strong 
subsidence modifying an elevated duct. Profile (c) typically occurs as the result of an 
evaporation duct over a body of water. Over land this type of M profile is common 
during the summer when nocturnal radiation over moist ground leads to a temperature 
inversion and a rapid decrease in moisture with height (Battan, 1973). Both profiles (b) 
and (c) are of direct concern to the mariner whose radar antennas are frequently found 


Within the limits of these ducts. 
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Fig. 2.1 Ducting occurrences for typical M profiles. 


Ii. DESCRIPTION OF SOME IREPS OUTPUT PRODUCTS 


Fig. 3.1 is a representation of one type of [REPS product. The “environmental 
data list” is used primarily for checking the numerical values of the input data. The 
IREPS computed output values of dew-point depression, dN/dz values, M units and a 
single word description of the refractive conditions give an overall assesment of the 


atmospheric refractive conditions at each input level. 
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Fig. 3.1 IREPS environmental data list for R/V Endeavor 
launch number I. 


The propagation conditions summary, reproduced in Fig. 3.2, is an EM system 
independent visual display and plain-language narrative assessment of expected 
refractive conditions. Vertical profiles of refractivity N as well as modified refractivity 


M are accompanied by a diagram showing the presence and vertical extent of any 


existing ducts. Evaporation duct height and surface wind speed information 1s given 
along with brief statements concerning anticipated performance of surface-to-surface, 
surface-to-air and air-to-air EM systems. 

Propagating radar waves experience constructive and destructive interference as 
direct path rays coincide with rays reflected from the earth’s surface. The resultant 
amplitude of the radar wave is a function of the phase differences among the 
intersecting waves. As a result, radar lobes (constructive interference) are created with 
shadow zones (destructive interference) appearing in between (Halliday and Resnick, 
1981). Knowledge of where these lobes and shadow zones occur can become an 
inumediate tactical advantage in both offensive and defensive operations. IREPS 
version 2.2 produces radar coverage diagrams which displavs the positions of the lobes 
and the shadow zones for any specified radar system and target at any desired 
probability of detection (POD). 

Fig. 3.3 is an example of the IREPS radar coverage display product which also 
can be very useful in determining an EM systems maximum range capability. It 
depicts a specified EM svstem’s area of coverage on a curved-earth, range-versus-height 
plot. Varying PODs are indicated by variations in the shading of the lobes. A numeric 
listing of some of the parameters used to generate the display, along with the location 


and time; date information of the profile, are included at the bottom of this product. 
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Fig. 3.2 IREPS propagation conditions summary for 
R/V Endeavor launch number I. 
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Fig. 3.3) IREPS radar coverage diagram for 
R/V Endeavor launch number I. 
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IV. DATA SETS AND DATA PROCESSING 


A. DATA SET SELECTION 

Two separate data sets were used for this thesis study. Both sets of radiosondes 
provided high resolution data. The CAPTEX data provided atmospheric information 
every 20 m throughout the profile, whereas the FASINEX atmospheric data were 
scparated by as much as 70 m at the surface. decreasing to 50 m near 500mb. A data 
set from a single station over land was chosen because the resulting analysis could be 
used to illustrate the evolution of atmospheric refraction conditions. The variability 
over land was expected to be very Tange uiiimeemmmacmee that over thesoceanas fic 
oceanic data were chosen to examine the notion that the MABL could be assumed to 
be homogeneous. Two ships in varying locations and separation distances were used 
to analyze the spatial vanabilty of the MABL. 

1. Cross Appalachian Tracer Experiment (CAPTEX) 

The high resolution overland balloon soundings used for this study were 
recorded as a part of the CAPTEX research program. It was conducted in response to 
the increased concern over regional and international aspects of air pollution and the 
recognized need for more reliable long range predictions of atmospheric pollutant 
transport and dispersion. The CAPTEX field program was designed to “provide data 
to test, evaluate and improve the transport modules of long range pollution models. 
CAPTEX should also provide insight into dispersion mechanisms. the role of vertical 
shear, mixing processes in and above the planetary boundary layer, the adequacy of 
standard meteorological observations for determining dispersion, and the role of the 
Appalachians and the Great Lakes in modifying synoptic and subsynoptic atmospheric 
processes’ (Ferber and Heffter, 1983). An extensive observation matrix, including 
supplementary (to the National Weather Service) rawinsonde launch sites and 
additional (up to four) launches each day, was designed to provide the data base 
necessary for forecasting and following the paths of released gaseous tracers. The 
CAPTEX data used in this thesis study were recorded at the Pennsylvania State 
University (PSU), which was one of the designated special rawinsonde sites. PSU was 
sct up to launch up to four soundings per day beginning on the tracer release days and 


for several days afterward. The special soundings were scheduled for 0200L, O800L, 
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1400L and 2000L (EDT). On the “non-experiment” days, two soundings ver day were 
iam iceditaesel pecs On O8YIL, 12 September 1983 (J.D. 255) and continues 
topougm Ul Oil sU October 1983 (J.D2 302), 

Prior analysis of these data for purposes not related to air pollution (Vaucher, 
1986) has clearly illustrated the significant variability in the atmospheric structure. 
Figs. 4.1, 4.2, and 4.5 show the time-height series of pressure altitude, temperature and 
Water vapor pressure, respectively. Marked variations in the temperature and 
especially the water vapor pressure as a function of time and height produced 
significant variations in the atmospheric refractivity, N, structure. Fig. 4.4 shows both 
optical and microwave mean refractivity (N units) profiles. The magnitudes decrease 
approximately exponentially with height, corresponding to what would be expected 
climatologicaily (Bean and Dutton, 1966). Fig. 4.5 illustrates the variabilitv of 
refractivity as a function of altitude: the standard deviation of N as a function of 
neight was plotted. Water vapor pressure is the dominant variable forcing the 
fluctuations of the microwave refractivity which are nearly 2.5 fold greater than those 
for optical refractivity. The relative maxima in standard deviation near the surface and 
at 2000 meters height are most likely due to changes in the height of the inversion layer 
and in the atmospheric refractivity in response to diurnal variations in radiative 
heating, nocturnal cocling, and the larger scale synoptic variability occurring during 
the period of the experiment. The relative minimum between the two can be related to 
the characteristic of the mixed layer which tends to keep conditions well-mixed 
resulting in smaller fluctations from the mean value of N than can be seen in cither the 
surface layer below or tne inversion layer above. 

2. Frontal Air-Sea Interaction Experiment (FASINEX) 

Previous air-sea interaction experiments, including JASIN (Pollard, 1978), 
have done little to determine what role horizontal variability might play in air-sea 
interaction. “However, the JASIN results were sufficient to point out the shortcomings 
of air-sea-interaction experiments that concentrate on observing the evolution and 
modification of boundary layer structure in only the vertical and to stimulate plans for 
experiments like FASINEX in which attention is focused on horizontal as vertical 
variability” (Stage and Weller, 1985). 

FASINEX was a “studv of the response of the upper ocean to atmospheric 
forcing in the vicinity of an oceanic front, the response of the lower atmosphere in that 


vicinitv to the oceanic front, and the two-way interaction between the ocean and the 
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atmosphere in the subtropical convergence zone southwest of Bermuda” (Stage and 
Weller, 1986). Although the entire experiment was conducted from the winter of 1985 
through the summer of 1986, the data used in this thesis study were obtained from 
radiosonde launches made by two oceanographic research vessels. The R/V Oceanus 
and the R V Endeavor made extensive meteorological observations throughout phase 2 
of the experiment during February and March 1986. The information gathered on 
board these two ships was used to investigate the horizontal homogeneity of the 
atmosphere. [he objective was to compare the differences in the atmosphere over each 
vessel while they remained, at times. on opposite sides of an oceanic front. 

In order to minimize contamination of the data set caused by temporal 
Variations in the atmosphere, only radiosonde launches that were made at similar times 
by each ship were compared. Of the 132 launches made by both ships, 13 pairs of 
radiosonde launches were within 37 minutes or less of each other. Table | gives the 
positions of each ship during the 13 pairs of launches, their separation distances, and 
the difference in radiosonde launch times. Figs. 4.6 and 4.7 display the position of 


each ship at the time of each launch. 


B. DATA PROCESSING 
1. CAPTEX Data 

Fig. 4.8 1s a flow digram of the CAPTEX data processing steps. The 
CAPTEX data were obtained from PSU (Thomson, 1986). In order to compare the 
differences in IREPS products on atmospheric soundings of differing resolution, the 
120 high resolution CAPTEX soundings were first degraded to resemble standard U.S. 
Navy operational soundings by selecting every fourth data point. This criterion 1s 
based on information from the Federal Meteorlogical Handbook No. 3, Radiosonde 
Observations (FMH #3). These two sets of CAPTEX data were then analyzed 
Separately to avoid prejudice in significant and mandatory level selection. After 
plotting the 120 temperature and dew-point temperature profiles of the degraded 
CAPTEX soundings on a standard skew T-log P chart, they were given to a qualified 
Aerographers Mate along with a copy of the pertinent pages of the Federal 
Meteorological Handbook No. 4, Radiosonde Code (FMH #4). AG2 Marty Schy then 
determined the significant and mandatory levels for each degraded sounding using the 
guidance in FMH #4. The temperature and dew-point profiles of the 120 high 
resolution CAPTEX soundings were plotted with the aid of GRAFSTAT, an 


interactive plotting routine avaliable on the IBM 3033 mainframe at the Naval 
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WevuiGe. | 


POSITIONAL INFORMATION FOR BOTH SHIPS 
AT PAE TIVES*Or eACH PAtiiter 


| 
LAUNCH SHIP POSITIONS SEPARATION DELTA 
NUMBER R/V OCEANUS R/V ENDEAVOR (KM) TIME (MIN) 
1 27.4N 69.6W 28.2N 70. 5W 125.6 O 
2 28.1N 70.3W 28.5N 70. 5W 50.0 26 
3 28.8N 70.6W 28.8N 70. 2W 37.6 15 
4 28.2N 70.1W 28.9N 70.0W 69.9 26 
|S 29.0N 69.9W 28.2N 70.1W 87.7 27 
5 28.2N 69.7W 27.9N 69. 8W 36.6 15 
7 28.7N 70.1W 27.2N 69. 8W NG 19 
8 28.2N 68.4W 28.4N 68. 4W Wl. @ 1 
9 28.6N 68.1W 28.5N 68.1W 9.0 26 
! 10 28.8N 68.1W 28.9N 67.9W 23.7 34 
1 27.3N 69.3W 28.8N 67. 6W 237.9 35 
12 27.1N 69.6W 28.6N 67.3W 276.8 3 
i 26.9N 70.OW 28.8N 67.5W Bee 4 15 


Postgraduate School. Significant and mandatory levels were determined by the author 
using the guidance of FMH #4. Pressure, temperature and dew-point temperature 
information for all significant and mandatory levels from both CAPTEX data sets were 
input into the IREPS version 2.2 refractivity code. Output information included, 
dM/dz values between the levels; ducting occurrences; and duct level, thickness and 
strength information. Statistical analvsis of these numeric values included average duct 
strength, average duct thickness. and standard deviation of the duct thickness for both 
the high resolution and degraded resolution profiles. The results of the statistical 
analysis of IREPS output appear in the results section, Chapter five, of this thesis. 
2. FASINEX Data 
Fig. 4.9 is a flow diagram of the steps involved in the FASINEX data 


processing. After choosing the 13 pairs of soundings within 37 minutes or less of each 
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Fig. 4.6 R/V Oceanus radiosonde launch positions. 





Fig. 4.7. R/V Endeavor radiosonde launch positions. 


oi 


other, the height of the sounding reaching the greater altitude was truncated so that 
each sounding pair would extend to the same altitude. These 26 launches were then 
prepared for the IREPS code by plotting, once again using GRAFSTAT, the 
temperature and dew point temperature profiles. Significant and mandatory levels were 
determined using the guidance in FMH #4. Pressure, temperature, and dew-point 
temperature information from each selected level, along with sea-surface temperature 
and surface wind speed information, was then entered into the IREPS code to produce 
the refractive structure of the atmosphere over each ship. Ship separation distance was 
computed using a simple great circle distance algorithm. Other statistical 
nianipulations were similar to those applied to the CAPTEX data. Additionally, the 
information was entered into a modified IREPS version 2.2 radar coverage diagram 
code to obtain radar propagation and coverage information. This modification did not 
change any of the mathematical manipulations used to derive the radar coverage 
diagrams; it merely involved the addition of programing statements which could 
determine the farthest extent of a lobe. For the purposes of this thesis study, the 
extent of a radar lobe was defined to be the maximum range and corresponding height 
reached by that lobe. With these two quantitative measurements of lobe 


characteristics, meaningful statistics and comparisons between lobes could be achieved. 
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Fig. 4.8 Flow diagram of CAPTEX data processing steps. 
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hig. 4.9 Flow diagram of the FASINEX data processing steps. 
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Phos CAPHREN RESULTS 

The results of comparing the high resolution radiosonde data with the degraded 
data are shown in Table 2. As one might expect, the high resolution data produced a 
far more detailed refractive structure, revealing almost five times as manv ducts. 
Because the high resolution data produced thin ducts in the profiles which were 
smeared out in the degraded data, the average thickness of the predicted ducts using 
the low resolution data was 31.8% larger than the average thickness of those from the 
high resolution data base. The standard deviation of the duct thickness was slightly 
larger as was the average duct strength for the high resolution soundings. This is due 
to the ability of the high resolution data to more accurately define the top and bottom 


of the ducts and to observe the greater M value differences between those heights. 





TABLE 2 


STATISTICAL RESULTS OF HIGH VS. DEGRADED RESOLUTION 
DATA 
HIGH DEGRADED 
| RESOLUTION DATA RESOLUTION DATA 
NUMBER OF DUCTS 68 14 
AVG. THICKNESS (M) GREE 152.79 
¢ DUCT THICKNESS (M) 57.58 50. 87 
AVG. DUCT STRENGTH 4.53 4.27 


(M UNITS) 





Another interesting result of comparing the degraded and the high resolution 
ducting layer results was that the high resolution sounding enabled resolution of 
multiple ducts in a single sounding. In some cases, the degraded sounding merely 
merged two slightly separated ducts into a single duct. This result is particularly 


important when considering that the minimum frequency that can be trapped by a duct 
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is a direct function of the thickness of the duct (Kerr, 1951). In most cases where a 
secondary duct was relatively shallow, the degraded resolution soundings missed the 
duct altogether. Fig. 5.1 1s a matrix of the number of ducts produced by the high 
resolution soundings compared to the number of ducts produced by the degraded 
resolution soundings for each launch. The strong bias toward the lower side of the 
diagonal is a direct indication of a high resolution sounding’s ability to better delineate 


the smaller scale refractive structure. 
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Fig. 5.1 Number of ducts discovered by each data set. 


Figs. 5.2 and 5.3 are histograms of duct heights for both the high resolution and 
degraded resolution CAPTEX data, respectively. This representation of the CAPTEX 
results once again shows that high resolution input data are required to observe the 
fine refractive structure of the atmosphere. Figs. 5.4 and 5.5 are a surface plot 
representation of duct evolution over time for both the high and low resolution data, 
respectively. Only the high resolution data are capable of realistically representing the 


ambient ducting conditions. 
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Fig. 5.2 Ilistogram of CAPTEX high resolution data duct heights. 
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Fig. 5.3 Histogram of CAPTEX low resolution data duct heights. 
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Visual inspection of the time variability of ducting features depicted in Fig. 5.4 
gives an indirect indication of how quickly the refractive environment can change. Not 
only are ducting features varying in height and extent over relatively short time 
periods, but they are also appearing and dissappearing. While MABL refractive 
variability may not be as dramatic as this continental data, IREPS users need to 
consider the timeliness of the input data, with respect to changing environmental 
conditions, before assuming atmospheric refractive conditions outputs are valid 


products capable of representing the present state of the atmosphere. 


B. FASINEX RESULTS 

Table 3 denotes the sea-surface temperature recorded at each ship during the 13 
pairs of launchings and, using an absolute difference in water temperature of at least 
one degree C, indicates whether the ships were on the same side or on opposite sides of 
an oceanic front. 

Using this information, Fig. 5.6 and Table 4 were derived. 

Distinct differences in MABL structure can be expected on opposite sides of an 
oceanic front (Stage and Weller, 1985). One would expect that the colder water side of 
the front should usually be associated with a more stable, less convective MABL and 
would thus allow higher evaporation ducts to exist. Conversely, the warmer water side 
Should enhance convective activity resulting in turbulent motion  supressing 
evaporation duct gradients. Fig. 5.6 shows a trend of the higher ducts existing on the 
colder side of the front supporting the above assumption; however, more data points 
are needed to precisely define this phenomenon. 

In an effort to ascertain if elevated duct occurrence is related to warm side/cold 
side variations, Table 4 provides some insight. The cold side elevated ducts were 
48.2% greater in mean strength and 51.2% greater in mean thickness than the warm 
side elevated ducts. Ducting occurences were 40% more prevalent on the cold side of 
the front. Again, since only six data points were used, the exact cause-effect 
relationships cannot be unambiguously determined; however, the results show a 
definite trend of larger, stronger and more frequent duct occurrences in colder water 
regimes. 

Table 5 shows the variability in elevated duct occurrences and evaporation duct 
heights observed by each ship. Fig. 5.7 shows that the differences in the evaporation 
duct heights observed by each ship are not a simple function of distance. The results 


of this figure indicates that one cannot assume similar atmospheric refractive 
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TABLE 3 


| 
OBSERVED SEA SURFACE TEMPERATURES DURING THE 
| FASINEX SHIP LAUNCHES 
| 
| 


LAUNCH SEA-SURFACE TEMPERATURES (°C) 
NUMBER R/V ENDEAVOR R/V OCEANUS 
1 Peo aCe) eee) 
| 2 Ui (ON 23.4 (W) 
| 3 Palas Dig 
4 Die 20 
5 23 OCW) Die aes 
6 2358 23.0 
7 2a 22.6 
8 papel 22.9 
9 Ving, (2 Dans 
| 10 poe 2 2G 
| iL DD (ome (2) 23.5 (W) 
| IL, 2 lle aes) 24.2 (W) 
i3 DDS 24.1 (WwW) 


(C) and (W) indicate cold and warm water areas respectively 





conditions based merely on short separation distances; conversely, great separation 
distances do not always imply large variations in refractive conditions. 

Using IREPS coverage diagrams, the variations in radar lobe heights and ranges 
between the pairs of launches were investigated. The most widely used radar in the 
U.S. Navy is the SPS-10 surface search radar found on board most surface combatants, 
amphibious ships and auxiliaries (Polmar, 1985). The SPS-10 radar parameters needed 
to generate the radar coverage diagrams are unclassified. Accordingly, it was the radar 
of choice for this trial. Comparing the differences in the maximum heights and ranges 
of the four radar lobes using a 50% POD on a one square meter fluctuating target, 
Figs. 5.8 through 5.15 again show differences in refractivity and radar propagation 


conditions of the atmosphere cannot be directly correlated simply with distance alone. 
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Fig. 5.6 Cold side versus warm side duct heights. 


[t is believed that only with knowledge of the synoptic-scale and mesoscale variability 
of the surrounding environment can meaningful interpretations of IREPS output be 
achieved. A real tactical disaster could result from blindly assuming homogeneous 
refractive conditions. Importantly, these results only depict refraction variations in a 
straight line path between two ships. Actual tactical situations involve 4-dimensional 
problems and the variations in atmospheric refractivity over the domain will be a 


function of the varying atmospheric conditions. 
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TABLE 4 


ELEVATED DUCT SUMMARY FOR WARM/COLD SIDE 
OCCURRENCES 
WITIT ELEVATION IN METERS AND STRENGTITIN M UNITS 


LAUNCH ELEVATED DUCTS 
NUMBER WARM SIDE COLD SIDE 
ELEVATION STRENGTH ELEVATION STRENGTH 
iei2. 65-1620. 52 22.24 
640.61- 810.28 6.81 
eos. 02-1714. 35 2.46 
2173.17-2243. 29 QO. 16 
1000. 09-1123. 6. 1120. 36-1444. 83 
1261. 05-1434. 
1269. 88-1381. Se 1190. 63-1590. 62 
63/./7/7= 836. : 1006. 46-1185. 62 
Zaso, 71-2939. 
Number of Ducts 
Average Thickness = 


Mean Strength 
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ELEVATED DUCT OCCURRENCES ANDEYArORAIION 
DLCT HEIGHTS: Ad BOTH SHIPS FORSEAGrigts UNCH 


LAUNCH ELEVATED DUCTS EVAPORATION DUCE BEICHE (uy 
NUMBER R/V ENDEAVOR R/V ENDEAVOR 
R/V OCEANUS R/V OCEANUS 
1 AL 0) ili eae 19.4 
ie Z 6) Zee oe 2 
3 0 6) 76 11.4 
4 0 6) ao A852 
> Z L Wats Zo 
6 al 6) 40.0 Za 
a 1 at Zoe 205 
8 0 Z cimas; 3 
9 1 1 1alseals: Pali s, 
10 0 1 ro, IES) 
ial 1 Z 10. 4 270 
1 ue 1 ACgS Zoes 
ES ul Z L055 WS. 3 
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Fig. 5.7. Differences in evaporation duct height as 
a function of ship separation distance. 
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Fig. 5.8 Differences in 50% POD lobe | heights for a SPS-10 radar. 
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lig. 5.10) Differences in 50% POD lobe 3 heights for a SPS-10 radar. 
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Fig. 5.11 Differences in 50% POD lobe 4 heights for a SPS-10 radar. 
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Fig. 5.12 Differences in 50% POD lobe 1 ranges for a SPS-10 radar. 
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Fig. 5.13 Differences in 50% POD lobe 2 ranges for a SPS-10 radar. 
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Pig. 5.14 Differences in 50% POD lobe 3 ranges for a SPS-10 radar. 
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Fig. 5.15 Differences in 50% POD lobe 4 ranges for a SPS-10 radar. 
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VI. POSSIBLE VARIABILITY IN LOCAL REFRACTIVITY 


A. SMALL SCALE REFRACTIVITY VARIABILITY 

Spaced-cavity refractometers supported by a captive ballon provide the ability to 
examine the small scale, atmospheric lavers of large variance in refractivity. Vertical 
Variations of as much as 8 N units in as little as one meter were discovered using this 
technique (Lane, 1964). The steepest gradients of refractivity were measured in 
inversion layers. Using aircraft in controlled descending and ascending flights, 
Fukushima and Inye (1964) found irregularities in refractivity exceeding 4 N units 
within 10 m in height. They also found middle-scale irregularities in the horizontal 
which they termed “refractivity clouds”. These refractivity clouds were characterized by 
a inedian intensity of 2 N units and a median scale of 600 m. Clearly, atmospheric 
“fine” structure 1s such that the concept of a well defined horizontally homogeneous 
layer is often an over simplification. 

Previous results of this study have concentrated on possible errors in atmospheric 
refractivity assessments caused by insufficient vertical resolution and by the 
assumptions of horizontal homogenity. The investigation into the effects of vertical 
resolution on duct detection suggests that the highest possible resolution is the most 
Gesirable. However, refraction variations in the ABL due to larger-scale turbulent 
eddies, as opposed to mean structure, could conceivably produce spurious ducts in an 
analysis if the data resolution is high enough. Organized structures, such as thermal 
plumes as first suggested by (Taylor, 1958), have been extensively investigated (e.g. 
Holland, 1968; Frisch and Businger, 1973; Khalsa and Businger, 1977; Lenschow and 
Stephens, 1980; Greenhut and Khalsa, 1982; Wilczak and Businger, 1983; and, Shaw 
and Businger, 1985). These investigations have shown significant variability in 
temperature and humidity, along with other atmospheric properties, between air inside 
and outside one of these features. Because the possibility of significant local variability 
in the atmospheric refractivity is real and could substantially affect IREPS analyses, 


some simple experiments using the CAPTEX and FASINEX data sets were conducted. 


B. CAPTEX POTENTIAL LOCAL VARIABILITY 
Table 6 is a numerical tabulation of microwave refractivity using both the mean 


and standard deviation values at selected levels derived from the 120 profiles taken 


50 


during the 47-day period of the CAPTEX experiment (Vaucher, 1986). In addition, 
two maximum deviation refractivity profiles were constructed by adding and 
subtracting half of the standard deviation values of refractivity from each of the mean 
values. A Monte-Carlo type experiment was designed such that six atmospheric 
profiles were constructed by randomly either using the value of mean N plus half of the 
standard deviation or the mean N minus half of the standard deviation for each level. 
i Oemicmre Milieu prolles \erewcmtercdmmto the IREPS code to produce their 
corresponding atmospheric refractivitv structure. Table 7 shows the results of some 
statistical manipulations performed upon the standard deviation derived profiles. 
Clearly, elevated ducting conditions are indicated quite frequently in the standard 
deviation-modified profiles whereas no ducts occurred from the use of the 


climatological mean profile. 


C. FASINEX POTENTIAL LOCAL VARIABILITY 

The previous experiment does not show the atmospheric variability associated 
with the intermittent larger scale-features. In order to investigate this possibility a 
second type of Monte-Carlo experiment was conceived. In this test, MABL properties 
reported in recent journal articles were applied to a typical FASINEX profile, which in 
turn, Was input into the IREPS refractivity program to determine possible variations in 
MABL refractivity structure that could be caused by the path of a radiosonde going in 
and out of these larger-scale features existing in the MABL. 

For all atmospheric temperatures, the second term and third terms of (2.3) can be 
combined so as to form (6.1) (Battan, 1973). This equation can be differentiated to 
determine the approximate effect of temperature and humidity variations on refractivity 


as shown in Equation 6.2. 


Ni= 77.6 P/T - 3.73 * 10° eT (6.1) 


AN = (6N,Gq)Aq + (6N/OT)AT (6.2) 


Shaw and Businger (1985) found that temperature variations in the open-ocean, 
extratropical MABL were insignificant (ie. 0.05 deg C). However, humidity 


fluctuations were of sufficient magnitude to affect MABL refractivity. Using Air Mass 
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TABLE 6 


CAPTEX MEAN AND STANDARD DEVIATION REFRACTIVITY 
Walle Jes 
FOR MICROWAVE PROPAGATION 


HEIGHT (M) MEAN N oN MEAN + oN/2 MEAN - oN/2 
0 318.0 14.9 325.5 311.5 
100 309.8 15.1 317.3 302.3 
200 305.1 14.3 ALY, 2 298.0 
500 294.3 127 300.6 288.0 
700 287.3 2a 293.4 JigieD 
800 Dosey, Ws gal 289.7 27a 
1000 275.5 oe 7 281.8 269.2 
1100 27250 12.8 278. 4 265.6 
1200 268. 2 13.3 274. 8 261.6 
1300 264.0 13.0 DIO), & 257.5 
1400 260.1 eee 266. 8 253.4 
1500 256.6 13.4 263.3 249.9 
Jere 248.3 ibs, 7 255.1 241.5 
1300 244.0 ete 250.5 237.5 
2000 236.8 125.5 243.0 230.6 
2400 2 eae 10.8 228.8 DME 
2700 214.50 9.72 219.36 209. 64 
2800 a ae 9. 43 215.96 206. 52 
3000 206. 41 9.22 Dimieg 201. 80 
3300 198. 42 8.08 202. 46 194. 38 
| 3500 193. 62 7. 88 197.56 189. 68 
3600 191.08 7.64 194. 90 187.26 
3800 185. 93 6. 80 189. 33 182.53 
4100 179. 26 6.20 182.36 176.16 
4300 162. 66 ay By 164. 85 160. 48 


S2 





Eee 


RES@GEVONS STATISTICS FROM CAPTEX STANDARD DEVIATION 
DE 2 D PROPILES 


Average number of ducts per profile Zig 
Standard deviation of duct thickness (m) ZO 
Average strength of ducts (M Units) ZO 


| 
| 
Average thickness of the ducts (m) IEA Frees: 
| 


Transformation Experiment (AMTEX) data, Lenschow and Stephens (1980) 
investigated the atmospheric structure of thermal versus non-thermal air in unusually 
convective regimes associated with cold air outbreaks over the Kuroshio current. Air 
inside these thermals can be identified by its upward vertical motion, humidity excess, 
and temperature excess (in the lower part of the MABL) when compared to air outside 
of a thermal. For refractivity calculations, increases in both temperature and humidity 
are compensating effects. Because we were interested in investigating a worst case 
scenario, We used information and equations from Lenschow and Stephens (1980) 
applied only to the humidity component of (6.2). Endeavor launch #1 was choosen as 
a typical FASINEN profile. The inversion layer height z, was found by inspection to 
be approximately 1700 meters. (6.3) describes the humidity excess inside a thermal. 


Sasa humidity scaling factor, and (6.4), 


Using an average value of mx of 0.095 g’m 
which compares the temperature excess inside ©, and the temperature excess outside 
©. a thermal, from Lenschow and Stephens (1980), (6.5) from (Fleagle and Businger, 
1980), ( 6.6) and (6.7), and a modified equation of state, (6.8), the possible variations in 
refractivity, AN, due to possible humidity variations at each level were computed. 
Since the inversion layer caps the MABL, the levels above the inversion would not be 
affected by variations in the MABL. Hence, the upper levels of the profile were not 
modified in this experiment. Once the magnitude of AN for cach layer was found, the 
layer was modified by  AN/2 to simulate possible variations in refractivity that could 
occur as a radiosonde passed into and out of thermals. As done previously with the 
CAPTEX data, six experimental profiles were generated by randomly adding or 


subtracting the value of AN/2 derived for each layer to the actual value of refractivity 
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for that layer. The resulting profiles were then entered into the IREPS 2.2 refraction 
code to see what effect, if any, would the intermittant larger-scale atmospheric features 
have on refractivity. Table 8 lists the resulting profiles which were subsequently used 
as input for IREPS. 


G,. = GES (7izaa se + 2(z/z,)°) ) (6.3) 
OR 0) 3,8), (6.4) 
LOS gicue lee a (6.5) 
Ride e7e7= 100 (6.6) 
q = .622 e/P (6.7) 
P = PRC ce eeiap (6.8) 


The only effect on the control profile was to modify the already existing trapping 
laver in thickness by plus or minus 17 meters, 5%, and its strength of 22.2 M units by 
plus or minus 2.1 M units, 9%. In two of the six cases, a weak surface trapping layer, 
dM/dz = -4 Km‘), appeared where previously a near trapping layer, dM/dz = 13.2 
Km), existed. The important result from this experiment is that temporary thermal- 
scale features in the MABL will have little effect on the mean refractive structure. 
Thus, spurious ducts from high-resolution measurements should not be a serious 


problem. 
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: | 
PELE & | 
MONTE-CARLO GENERATED N PROFILES FROM A TYPICAL 
FASINEX SOUNDING 
: 
| 


Pe cola oi notezet Bot 3 TEST 4 TEST 5 TEST 6S 
oS. 0 343.0 343.0 343.0 343.0 343.0 343.0 
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3 S253 O26 Omer omommc Zo wom SZS.6 323.6 321.4 

“ SAN, Ce oreo S22 622.48 S228 S22. 1 
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oe 274. 2 if Crome omc 62/6. 3 Zien 27241 


* Level 6 represents the top of the MABL 





This same type of experiment was performed with the CAPTEX data. The 
temperature variation instead of the humidity variation term of Equation 6.2 was used 
for two reasons. First, humidity variations over land are generally smaller than 
hunudity variations over water. Secondly, we wanted to compare the relative effects of 
humidity and temperature variations on refractivity. After applying the same 
procedure to a typical CAPTEX profile as was done previously with R/V Endeavor 
launch number one, we found virtually no effect on the control N profile from the 
surface to the ABL. Lane (1964) confirms this result by showing that fluctuations in 
humidity mixing-ratio have approximately a seven times greater effect on refractive 


index than do temperature fluctuations. 


a0 


Vil. SUMMARY AND CONCLUSIONS 


IREPS version 2.2 is the newest refractive effects prediction system to be used by 
the Fleet. While it will provide important tactical information to operational 
commanders, the system will be limited in its ability to reflect the actual atmospheric 
refractive conditions by the quality and timeliness of the data which is currently used 
as input. By comparing the [REPS version 2.2 computed atmospheric refractivity 
conditions for a high resolution and degraded resolution sounding, we showed the 
importance of using high resolution data. The high resolution data were better able to 
define the smaller scale refractive structure of the atmosphere and, therefore, allowed 
IREPS to realistically portray the ambient ducting conditions. 

The high resolution CAPTEX data also provided insight into how quickly 
atmospheric refractivity conditions can change. By observing the variabiltv in ducting 
conditions for 120 soundings taken at irregular time intervals over a 47-day period, we 
found that refractive conditions are in constant dynamic evolution on diurnal as well as 
synoptic time scales. While refractive variability may not be as pronounced over water 
as it 1S Over a continent, the implication that timely data must be used to achieve 
accurate refractivitv condition predictions still should be considered by all IREPS users. 

Although the IREPS version 2.2 user’s manual implies horizontal inhomogeneity 
of the atmosphere should not be a serious refractivity consideration 85% of the time, 
the results of the investigations using the FASINEX data showed horizontal 
inhomogeneity does indeed cause significant variability in radar lobe coverage at least 
50% of the time. This variability is not a function of ship separation distance, but 
appears to be driven by atmospheric conditions surrounding both platforms. 

The last section of this thesis investigated the possibility of intermittently 
existing, relatively large-scale, atmospheric features creating spurious ducts because of 
the use of high resolution data. By applying worst case scenarios to our data base and 
entering modified sounding profiles into the IREPS version 2.2 refractivity code, we 
could compare the resulting output with the output of the original sounding. In both 
cases, with CAPTEX and FASINEX data, there was no significant modification of the 
previously eXisting refractive conditions. Accordingly, we see no disadvantages to 


using high resolution data. Quite to the contrary, we believe that only by the use of 
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umely, high resolution sounding data can an accurate assessment of atmospheric 
refractive conditions be derived by IREPS version 2.2. 
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